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ABSTRACT 



We have discovered an extended Lyman-a plume associated with a compact 
source at redshift z 5.4 in slitless spectroscopic data from the Grism ACS 
Program for Extragalactic Science (GRAPES) project. The spatial extent of 
the emission is about 6 x 1.5 kpc (1" x 0.25"). Combining our grism data and 
the broadband images from the Hubble Ultra Deep Field (UDF) images, we 
find a Lyman-a line flux of ~ 2.2 x 10~^^ erg cm~^ s~^ and surface brightness 
~ 7 X 10~^^ erg cm~^ s~^/n". The UDF images show diffuse continuum emission 
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associated with the Lyman-a plume (hereafter UDF 5225), with three embedded 
knots. The morphology of UDF 5225 is highly suggestive of a galaxy in assembly. 
It is moreover possible that the prominent Lyman-a emission from this object is 
due to an active nucleus, and that we are seeing the simultaneous growth through 
accretion of a galaxy and its central black hole. FoUowup observations at higher 
spectral resolution could test this hypothesis. 

Subject headings: galaxies: high redshift — galaxies: formation — galaxies: 
interactions — galaxies: starburst — galaxies: individual (UDF 5225) 



1. Introduction 

The Grism ACS Program for Extragalactic Science (GRAPES) project^ is a slitless 
spectroscopic survey that exploits the potential of the G800L grism on the Hubble Space 
Telescope's Advanced Camera for Surveys (ACS) to achieve the most sensitive unbiased 
spectroscopy yet. GRAPES is targeted in the Hubble Ultra Deep Field (UDF) region, to 
complement the UDF direct images, which are in turn the deepest optical imaging to date 
(Beckwith et al. 2004). The GRAPES survey, and in particular our data analysis methods, 
are described in more detail by Pirzkal et al. (2004). 

One of the primary scientific goals of GRAPES is to study the luminosity function of 
Lyman break galaxies (LBGs) using spectroscopically confirmed samples at unprecedented 
sensitivity, and thereby to constrain the faint end luminosity function slope. We have begun 
this effort through a targeted look at photometrically selected Lyman break candidates, 
using both i-dropout galaxies from the UDF (Malhotra et al 2004) and F-dropout galaxies 
which we have identified (Rhoads et al. 2004) following the selection criteria outlined by 
Giavalisco et al. (2004b) in both the the vl.O GOODS survey data (Giavalisco ct al. 2004a) 
and the UDF. (We refer to ACS and NICMOS filters by names of roughly corresponding 
ground-based filters: F435W ^ B; F606W ^ V; F775W -> i; F850LP ^ z; FllOW 
J; and F160W — > H.) Most of the confirmed Lyman break objects are spatially compact, 
with sizes (< 0.5") and morphologies typical for the LBG population (e.g., Ferguson et al. 
2004). In this Letter, we describe the most prominent exception to this pattern we have 
encountered to date, a F-dropout object (designated UDF 5225) that is exceptional in its 
size, morphology, and spectroscopic properties. 



^Accepted for publication in The Astrophysical Journal, with tentative publication date March 2005 
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Throughout this paper we use the current concordance cosmology {Hq — 71 km s~^ Mpc~^ 
= 0.27, ^totai — 1; see Spergel et al. 2003). Magnitudes are given on the AB 
magnitude system, so that magnitude zero corresponds to a flux density /j, = 3.6kJy = 
3.6 X 10^^° erg cm~^ s~^Hz~^. We present the observations in section 2, compare them to 
various physical models for UDF 5225 in section 3, and summarize our conclusions in sec- 
tion 4. 



2. Observational Properties of UDF 5225 

UDF 5225 has two morphological components (see fig. 1): A compact "core," which is 
barely resolved in the ACS data (with full width at half maximum ^ 0.12" in the UDF z 
filter image, which has a point spread function [PSF] with FWHM ^ 0.10") and an extended 
"plume" with a size ~ 1.0" x 0.3". There are three distinct condensations or "knots" within 
the plume. The optical and near-IR colors of the object (discussed below) identify it as 
a Lyman break object, with an intrinsically blue spectrum suppressed by the inter galactic 
medium at V band and bluer wavelengths. 

The size is unusual relative to 2; ~ 5 Lyman break galaxies selected photometrically 
from the GOODS data (Ferguson et al. 2004), which show a broad peak between 0.1" and 
0.5" (comparable to the minor axis size of UDF 5225) and no galaxies as large as 1" (the 
major axis size of UDF 5225). UDF 5225 itself is not in the Ferguson et al. sample, being 
slightly fainter than their flux limit. 

The 2D ACS grism spectra of UDF 5225 are shown in fig. 2. They detect UDF 5225 
significantly in each of the five epochs analyzed (see Pirzkal et al. 2004, and see also Riess 
et al. 2004 for more detail on epoch 0). Our strategy of using many roll angles results in a 
clean separation of the core and plume spectra for epochs 0-2 [PA = 117°, 126°, and 134°, 
where "PA" refers to the position angle of the Hubble Space Telescope's V3 axis), while the 
spectra from the two components are superposed in epochs 3 and 4 (PA = 217° and 231°). 
The spectrum from the northwestern tip of the plume is contaminated by the spectrum of 
an unrelated, brighter object in the PA = 134° data, but otherwise the UDF 5225 spectra 
are free of significant overlap. 

Where the plume's spectrum can be examined independently of other sources, including 
the core (i.e. epochs 0, 1, and in part 2), it is dominated by a single strong emission line at 
~ 7800A. Because this line falls in the i filter bandpass, and the plume is detected in both 
i and z images, we know that there must also be weak continuum emission from the plume 
on the red side of the hne. The core shows both a break and a hne at the same wavelength. 
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When both core and plume component spectra are superposed (epochs 3-4) , their combined 
hnc and continuum flux results in a stronger spectroscopic detection. We identify the line 
and break with Lyman-a, based on their wavelength coincidence in the spectrum of the core 
and based on the i?-band nondetection and very weak V^-band flux of the source. This then 
implies a redshift z ~ 5.42, with an estimated uncertainty bz ^ 0.07. The object is near 
the upper end of the redshift range for V dropouts (and approaches the redshift range of i 
dropouts) . 

We measured the broadband optical magnitudes of the core and plume components 
using the UDF images. We deflned apertures to match the morphology of the core and plume 
components. These apertures follow isophotes in a version of the UDF i band image smoothed 
with a 0.14" FWHM Gaussian kernel, except at the boundary between the two components. 
There is no deep minimum in surface brightness separating the core from the plume. We find 
for the core i = 27.94 ± 0.017, i-z^ 0.51 ± 0.025, andV-i^ 2.3 ± 0.12 mag, while for the 
plume we obtain i = 27.39 ± 0.023 mag, i — z = 0.45 ± 0.034 mag, and V — i = 3.3 ± 0.45 mag. 
Comparing our i-band flux with the flux in the UDF catalog from STScI, we flnd that the 
masks contain about 55% of the overall flux of the source, with the remaining flux in lower 
surface brightness regions that we did not use for our color measurements. We did not 
attempt to measure the colors of the knots within the plume individually, but inspection of 
figure 1 shows that they are consistent with the overall colors of the plume. Comparing the 
colors of the main components, we see that the core is slightly redder ini—z and considerably 
bluer inV—i, but that the significance of both statements is low: About 1.5(7 for i — z and 3a 
iorV — i. (Note, we calculated the significance of the V — i color difference directly from fiux 
ratios rather than magnitudes.) The V — i colors are consistent with transmission through 
the Lyman-a forest of photons emitted at 912A < Ajcst < 1215A. For a flat {f^ ^ constant) 
continuum, we would expect the Lyman-a forest to attenuate the V band flux by a factor of 
14 for a source at redshift z = 5.4, based on the formalism of Madau (1995). 

To eliminate any possibility that UDF 5225 is simply a faint, intrinsically red galaxy 
at a lower redshift, we also examined the near-infrared colors of the object, using NICMOS 
images of the UDF (PI R. Thompson). We flnd a blue color {z — J = 0.04 and J — 
H = —0.27). The NICMOS measurements were performed using an elliptical aperture large 
enough to encompass the entire object (size «i 1.0" x 0.25"). We did not attempt separate 
IR measurements of core and plume colors, given the larger pixel size, point spread function 
size, and lower sensitivity of the UDF NICMOS data. The observed NIR colors support 
the identification of UDF 5225 as a Lyman break object with an intrinsically blue spectrum 
truncated below Lyman-a by the intergalactic medium. 

Because UDF 5225 is both faint and extended, it is near the practical detection limit 



- 5 - 



of the GRAPES data. We therefore combine the UDF imaging data, the NICMOS UDF 
observations, and the GRAPES redshift to estimate the hne flux and equivalent width by 
fitting the spectral energy distribution (SED). We model the intrinsic spectrum as a power 
law continuum plus an unresolved Lyman-« emission line, modifying both by the Lyman-a 
forest transmission calculated under the model by Madau (1995). The z — J and J — H 
colors are unaffected by Lyman-a emission and Lyman-a forest absorption, so we use them 
to constrain the intrinsic spectral slope to a ^ —2.4 ± 0.3, where fx oc A". With the 
slope fixed, the i — z and V — i colors are determined by the line flux and the redshift 
(which determines how strongly Lyman-a forest absorption affects broad band fluxes). For 
z — 5.42 ± 0.07, we find a rest frame equivalent width of 70 ± 30A and an observed line 
flux of (2.2 ± 0.8) X 10~^^ergcm~^s~^. The range in line flux and equivalent width is 
primarily determined by the range of acceptable redshifts. Changing the continuum slope 
within its plausible range has a rather smaller effect on the line flux needed to match the 
i — z color. The observed line flux corresponds to an approximate surface brightness of 
7 X lO-^^crgcm-2 s^VD". 

To test whether the Lyman-a emission in UDF 5225 could be powered by an active 
galactic nucleus (AGN), we examined the UDF direct images from multiple epochs for vari- 
ability. We stacked the UDF i- and 2;-band data into eight epochs of approximately equal 
exposure time. We also stacked archival z-band imaging from HST program 9352, which is 
shallower than 1/8 of the UDF data but gives a longer time basehne (~ lyr). Subtracting 
the mean of the UDF stacks from each individual epoch shows no significant residuals at the 
location of UDF 5225, leading us to conclude that the source is not significantly variable in 
the UDF data, with an upper limit to flux variations of < 10%. Unfortunately, at z ^ 5.4, 
any Civ line lies beyond the wavelength coverage of the grism (as do all redder AGN lines), 
while the Nv line is not separable from the Lyman-a line at the resolution of the grism. 
UDF 5225 is not detected in X-rays, based on the 1.0 Msec Chandra Deep Field South cata- 
log (Giacconi et al. 2002; Alexander et al 2003), though given its high redshift, this does not 
strongly exclude an AGN: Only one AGN at, z > 5 has been discovered in the two Chandra 
deep fields so far (Barger et al. 2002), and an X-ray detection would only be expected for 
a rest frame hard-band luminosity > 1.8 x lO^^ergs"^. Higher resolution optical spectra of 
UDF 5225 could convincingly determine whether or not it harbors an AGN by measuring 
the velocity width of the Lyman-a emission. 
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Fig. 1. — Color composite image of UDF 5225 from the Hubble Ultra Deep Field data. The 
elongated red source in the center is UDF 5225. We have labeled the "core" component at 
the bottom and the "plume" extending towards the top of the image. We also label the three 
knots in the plume as "Kl," "K2," and "K3.". The length of the plume is nearly 1". Blue 
colored objects nearby are unrelated foreground sources. Red represents the z filter, Green 
an average of V and i, and blue the B filter. Based on the color composite UDF image by 
Z. G. Levay. 
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3. Models for UDF 5225 

What is the nature of the Lyman-a plume in UDF 5225? The hne luminosity and the 
equivalent width are both broadly similar to those seen in narrowband-selected Lyman-a 
samples (e.g., Rhoads et al. 2000, Malhotra & Rhoads 2002, Ouchi et al. 2003, Hu et al. 
2004). The morphology is reminiscent of radio galaxies (e.g. Windhorst et al. 1998), quasar 
jets, the "Lyman-a blobs" that have been observed at 2; ~ 3 (e.g., Steidel et al. 2000), 
and some other high redshift galaxies (Pascarelle et al 1996, Bunker et al. 2000, Keel et al. 
2002). 

We can consider several possibilities: A recombination nebula powered by in situ star 
formation, possibly triggered by a galactic merger; light from the core component scattered 
by either electrons or dust; or a a recombination nebula powered by the core. 

In situ star formation: Star formation would provide a local source of ionizing photons 
within the plume. It also produces lower-energy photons that would dominate the rest-UV 
continuum redward of Lyman-a. The resulting continuum slope is the bluest among the 
models considered here. The rest frame Lyman-a equivalent width for star formation should 
be < 240A based on stellar population models at Solar metallicity (Chariot & Fall 1993). 
While observations of Lyman-a galaxies at z ^ 4.5 often show larger values (Malhotra & 
Rhoads 2002), the presence of stars in situ will always reduce the equivalent width relative to 
recombination models where the ionizing photon source is distant. Star formation is consis- 
tent with essentially any morphology. A conversion factor of 1Mq/jt — 10^^ ergs"-*^ is widely 
used for high redshift Lyman-a emission, and would imply star formation at > IMq/jy in the 
plume, subject to the standard (but untested) assumptions that the star formation follows a 
Kennicutt (1983) initial mass function (IMF) and that case B recombination is valid. If we 
base our estimate instead on the rest frame UV continuum emission at 1425A, as measured 
by the z band image, we find a star formation rate of ~ 4M0/yr for the plume (plus an- 
other 2.5M0/yr in the core, assuming the core light is not dominated by an active nucleus). 
The conversion factors for both continuum and line light are substantially uncertain, due 
to assumptions about the IMF and about the effects of gas and dust on Lyman-a radiative 
transfer. Thus, this constitutes remarkably good agreement. 

Scattered Light: If UV radiation from the core encounters a sufficiently dense scatter- 
ing medium, a detectable scattering cone could be produced. This is the least plausible 
explanation for the UDF 5225 plume, because the core luminosity would have to be very 
large to power electron scattering, while dust scattering would be unlikely to produce a large 
Lyman-CK equivalent width and a blue continuum color. 
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We estimate the mass of scattering material required under the approximation that the 
plume is a cone with length 6kpc and base diameter ~ 2kpc. The corresponding volume 
(in physical, not comoving, units) is 2 x lO^^cm^. We estimate the mass for a scattering 
optical depth Tg from the core to the end of the cone. Then for electron scattering (with 
(Tt = 1.6 X lO^^^cm^) we find a number density of ~ 25re and a corresponding mass of 
~ 5 X lO^TeM© in the cone. This gas would emit copious Lyman-ct radiation. Indeed, to 
avoid producing more than the observed Lyman-a luminosity would require that < 0.025. 
Clumping of the scattering gas would be required to reproduce the observed knots and the 
30° bend in the plume. This clumping would further enhance recombinations and further 
reduce the maximum Tg for electron scattering to dominate the plume emission. 

Such low optical depths imply that electron scattering is a very inefficient way of produc- 
ing the Lyman-CK nebulosity, because scattered light will be suppressed by a factor t^- The 
scattering region must then be illuminated at a much higher intensity than one would naively 
infer from the observed core flux. Consider a simple toy model where the central source has 
two emission components, one with a ~ 30° opening angle that powers the observed plume, 
and one isotropic that powers the observed core and that we allow to be attenuated by optical 
depth Tabs of absorption. To reproduce the observed plume to core flux ratio, the collimated 
component would then need to contain a fraction ~ 1/ [1 + 0.6re exp(Tab5)], which becomes 
~ 1 — 0.6re > 98% in the hmit where Tabs ^ 1- This would correspond to a total source 
UV luminosity vL^, > 3 x lO^^ergs"^ (measured in the z-band, which is Arest ~ 1410A) i.e., 
7 X IO^^Lq or an angle- averaged absolute AB magnitude of —23.7. This luminosity would 
increase if either the assumed or degree of coUimation were reduced. Such a model would 
likely require an AGN in the core, because star light cannot be tightly collimated, which 
would increase the luminosity requirement by another factor of > 10 while requiring also 
'Tabs ^ 2 on the line of sight to the core. The mass of ionized hydrogen required in an electron 
scattering scenario would be modest, < lO^M©. 

Dust scattering in an ionized medium could work, as long as the optical depth is suitably 
small so that (a) there is no significant reddening of the scattered continuum light, and (b) 
Ly-a radiation is not selectively absorbed relatively to the continuum. Condition (a) requires 
T < 0.1 at Arest ~ 1300A, and similarly condition (b) requires r < 0.05 at Arest ~ 1216A 
to avoid attenuation of Lyman-a by more than a factor of 2 (Panagia and Ranieri 1973a, 
b). The associated total mass of gas and dust would be Mpiume ~ lO^Tdust (1216A)Mq 
assuming a dust-to-gas mass ratio of O.Ol. Inserting Tdust < O.l gives -Mpiume lO^M©. The 
requirement on the core luminosity can be derived exactly as for the electron scattering case, 
but now with larger scattering opacity, so the final constraint becomes vL^, > 3 x 10^^ erg s~^. 

An additional concern with scattering models is that the intensity of scattered light 
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should decrease away from the AGN (falling off as in the simplest case, though projection 
effects and anisotropic quasar emission could modify this) . Yet the surface brightness of the 
plume is constant to within factors of two over a factor of ~ 8 in distance from the core, 
which would require the column density to increase with distance as ~ r^. 

Recombination powered by core light: Objects like "Lyman-a blobs" (e.g., Stcidcl ct 
al. 2000) and some radio galaxies have extensive Lyman-a nebulae (up to lOOkpc in size) 
that are likely powered by AGN emission. We now consider whether UDF 5225 could be a 
physically similar object on a smaller scale. 

Such a model is similar to the electron scattering model but with higher densities, where 
recombinations become more important than electron scattering. In this case the mass of 
ionized hydrogen involved would be > IO^Mq. The core luminosity could be much lower 
in this case: For an ionization bounded nebula, it need only be about twice the measured 
Lyman-a luminosity, i.e., > 2 x 10*^^ergs~^ in ionizing radiation, although it could be much 
larger if the radiation is largely isotropic and the morphology of the plume is set by where 
there is substantial gas. 

However, in this case, the continuum emission from the plume would be pure nebular 
emission, with components from the two-photon process, bound-free emission, and free-free 
emission. The expected equivalent width of the line would then be > lOOOA (rest frame) and 
the z band continuum should be much weaker than we see. While dust attenuation of the 
Lyman-CK might help, the model would have to be fairly contrived to simultaneously fit the 
Lyman-a flux, continuum flux, equivalent width, and color in the plume, and we therefore 
disfavor this model also. The only countervailing argument comes from the V band flux, 
which should be absent in a two-photon continuum and is indeed rather weaker in the plume 
than the core. However none of the V band detections arc strong, so it is not clear if the 
difference in V — i and V — z colors between the two components is significant. 

A Merger Scenario: Consider now the possibility that UDF 5225 is a major merger in 
progress at z ~ 5.4. The nuclei of the two interacting galaxies would then presumably be 
the "core" component we have discussed plus the brightest "knot" in the plume, which lies 
0.63" away. Quantitative morphological tests based on the asymmetry parameter A (e.g., 
Gonselice et al. 2003) are consistent with this scenario, based on an analysis of GRAPES 
object morphologies now in progress (Pirzkal et al. 2004b). If we measure asymmetry 
relative to the center of light for UDF 5225 (i.e., almost half way from the core to the end 
of the plume), we find A{z) = 0.34 ± 0.09 and A{{} = 0.16 ± 0.09 for the z and i filters. If 
we instead measure A with the center placed in the core component, we get much higher 
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numbers, A{z) — 0.45 ± 0.13 and A{i) — 0.43 ± 0.13. Mergers are typically found to have 
A > 0.35 (Conselice et al. 2003). Further discussion of these parameters for the full i- and 
V^-drop GRAPES samples will be presented in Pirzkal et al (2004b). 

The projected separation of the core and knot, ~ 4kpc, would then imply a crossing 
time of order 40t'iooMyr, where vwo ~ 1 is the relative velocity of the two components in 
units of 100 km s~^. Multiplying this by the star formation rates inferred in section 3 implies 
formation of some few x IO^Mq of stars in the course of the interaction. 

The star formation rate per unit area, based conservatively on the UV-derived SFR for 
the plume, is ~ O.4M0yr~^kpc~^. Comparing to the global Schmidt law for star formation, 
we would infer a gas mass surface density of ~ SOOA/qPC"^, for a total gas mass of order 
7 X IO^Mq. This number implies a gas consumption time scale of ~ Wt^yn ~ 600Myr, 
though if star formation is proceeding at an atypically high rate (driven by interaction), the 
gas reservoir could be smaller. 

4. Discussion 

We have examined several possible scenarios for the observed properties of the galaxy 
UDF 5225, a very faint and high redshift object with a core-plume morphology and prominent 
Lyman-a emission. We conclude that the Lyman-a emission is most likely powered by in 
situ star formation throughout the object. Present evidence neither requires nor rules out 
the presence of an AGN in the core component. FoUowup spectroscopy using higher spectral 
resolution and/or coverage into the near infrared would provide new information that could 
settle the AGN question. Tidally triggered star formation in a merging galaxy pair seems to 
describe the galaxy well. As such, it may be a particularly spectacular example of the broad 
class of star-forming Lyman break galaxies that dominate the galaxy population observed in 
the Hubble Ultra Deep Field and other very sensitive, high redshift galaxy surveys. 

We thank the STScI Director and the UDF team for their hard work in designing and 
executing the UDF experiment. We thank Zoltan G. Levay for producing the color composite 
shown in figure 1. This work has been supported under grant number HST-GO-09793. C.G. 
acknowledges support from NSF-AST-0137927. L.A.M. acknowledges support by NASA 
through contract number 1224666 issued by the Jet Propulsion Laboratory, California Insti- 
tute of Technology under NASA contract 1407. 
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Fig. 2.— The GOODS survey vl.O direct image of UDF 5225 in the z' (F850LP) filter (left 
panels), together with the 2-D GRAPES spectra at each of five roll angles (right panels). 
From top to bottom, the angles shown are PA 117, 126, 134, 217, and 231°. The direct image 
is the same in all rows but is rotated to match the roll angle of each 2D grism spectrum 
cutout. All images have been smoothed with a Gaussian filter matched to the angular size 
of the plume component (0.3" FWHM) for maximum clarity. The upper edge of the grism 
spectrum in the third PA shows contamination by an unrelated source elsewhere in the image. 
The Lyman-a emission from both plume and core is most clearly visible in the second and 
third panels; the first panel is at a similar roll angle but is shallower. In the fourth and 
fifth panels, the line emission from the plume falls atop continuum emission from the core, 
resulting in a higher surface brightness in the dispersed image. 



